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Abstract. Based on the path integral approach the partition function of a many body system with separable 
two body interaction is calculated in the sense of a semiclassical approximation. The commonly used 
Gaussian type of approximation, known as the perturbed static path approximation (PSPA), breaks down 
near a crossover temperature due to instabilities of the classical mean field solution. It is shown how the 
PSPA is systematically improved within the crossover region by taking into account large non-Gaussian 
fluctuations and an approximation applicable down to very low temperatures is carried out. These findings 
are tested against exact results for the archetypical cases of a particle moving in a one dimensional double 
well and the exactly solvable Lipkin-Meshkov-Glick model. The extensions should have applications in 
finite systems at low temperatures as in nuclear physics and mesoscopic systems, e.g. for gap fluctuations 
in nanoscale superconducting devices previously studied within a PSPA type of approximation. 



PACS. 5.30.-d quantum statistical mechanics - 24.60. 
thermodynamic properties - 21 TO. Ma level density 



k statistical theory and fluctuations - 74.25.Bt 



1 Introduction 

Thermodynamic properties of interacting many body sys- 
tems are of fundamental interest for all kinds of condensed 
matter. A particular challenge has been the temperature 
range where strong quantum effects render simple mean 
field theories insufficient. For finite systems with separa- 
ble two-body interaction much can be gained by relying on 
approximate techniques to evaluate the partition function. 
Typical examples include certain mesoscopic systems, e.g. 
ultrasmall superconducting metallic grains pi on which 
substantial research has focused recently. 

A very elegant approach to approximate the partition 
function is provided by the path integral formalism. There, 
the usual trick |||J i s to turn the two-body interaction via 
a Hubbard-Stratonovich transformation into terms con- 
taining only one-body operators and an auxiliary field. 
While the static part of this field determines the static 
mean field result, many body quantum fluctuations are 
encoded in its dynamical part. The latter one can sys- 
tematically be accounted for in the sense of a semiclas- 
sical approximation. Accordingly, a static approximation 
has been developed originally in order to study finite size 
effects in small superconductors j|. The static path ap- 
proximation (SPA) has also been used in nuclear physics 
to calculate thermodynamic properties and level densi- 
ties of hot nuclei Later it has been extended by the 
SPA+RPA the perturbed static path approximation 
(PSPA, the name which we take over in the following) 0, 



or the correlated static path approximation (CSPA) ||£§. 
These approaches account for small Gaussian fluctuations 
around the classical mean field and are thus applicable in 
a temperature range where quantum effects are no longer 
negligible. This way, from the approximate partition func- 
tion a variety of thermodynamic quantities have been cal- 
culated in the literature: the free and internal energy, the 
specific heat and the level density of the system |p|,[7|,[ic|, 

f. In addition thermal expection values of observables 
and strength functions MM have been deduced. While 
many of these calculations have been carried out in ex- 
actly solvable models ||@,|],[|] also reasonable agreement 
with experimental data was shown for the level density 



of finite nuclei 
put forward in 



lC,pT[. Further, the formalism has been 
11] to take into account the coupling to 
the continuum and in [|l2| an extension of the PSPA to 
describe dissipative decay out of metastable states has 
been developed in a self-consistent fashion. Very recently, 
the improvements of the SPA have been applied to small 
strongly correlated condensed matter systems. For exam- 
ple, using the CSPA odd-even effects in small superfluid 
systems |14 and motivated by new experiments on ultra- 
small superconducting metallic grains jlj gap fluctuations 
and pairing effects in finite size superconductors have been 
examined fig ). 

As noted above, these conventional improvements of 
the SPA take into account quantum effects on the RPA 
level only, i.e. Gaussian fluctuations, and thus, are plagued 
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by the problem that they break down at a certain tem- 
perature To- Physically at this temperature the classical 
mean field solution becomes unstable in functional space 
and large fluctuations render the Gaussian approxima- 
tion insufficient. This has some features in common jl3], 
|i2| ] with the change of stability at a crossover tempera- 
ture occuring in the description of dissipative decay e.g. 
in Josephson junctions |]l6| but is different in detail. By 
thoroughly analysing the appearence of instabilities for 
the many body problem at lower temperatures, we go be- 
yond the PSPA to achieve a smooth behavior around To. 
For temperatures sufficiently below To the partition func- 
tion is dominated by field configurations away from the 
classical mean field around which a Gaussian approxima- 
tion is again valid. In this paper we lay out the general 
theory and discuss its validity in comparison to exactly 
solvable models before it can be applied to realistic sys- 
tems. 

We start with a brief outline of the main results of the 
conventional PSPA and introduce the notation used in the 
sequel (Sec. ||). The extension of the PSPA around the 



number (A) fixed on average. In principle, one should work 
with truly fixed particle numbers, i.e. within the canoni- 
cal ensemble. However, as we are mainly interested in the 
dependence of system properties on excitation energy or 
temperature, it is more convenient to exploit (|^). A very 
elegant way to do this is to represent Z{(3) as a func- 
tional integral in imaginary time 0], which also allows to 
systematically include fluctuations around the mean field. 
Given the form of the Hamiltonian in (0) the mean field 
approximation starts with a Hubbard-Stratonovich trans- 
formation pO of the imaginary time path integral corre- 
sponding to (3). Accordingly, the product FF is split by 
introducing an auxiliary path q(r) as a collective variable. 
Since this procedure is well-known (see e.g. 0,§,||) we 
simply state here the basic results which will then serve 
as the starting point for our analysis ([]). 
After introducing the Fourier expansion 

q(r) = qo + ^2 q r exp(ii/ r r) , q_ r = q* , (3) 

r#0 



critical temperature T is developed in Sec. O which in with the Matsubara frequencies 



Sec. 2.3 allows for a treatment of the temperature range 
far below To. In Sec. || our results are first applied to 
a one-dimensional double well, before in Sec. |] we turn 
to the thermodynamic properties of the Lipkin-Meshkov- 
Glick model g7j. 



2 Partition function of an interacting many 
body system 

In this paper we want to approximate the partition func- 
tion of a system described by a Hamiltonian of the follow- 
ing structure 



2ir 2ir 

v r =-r= T rT , , = ±1, ±2, ±3 



(4) 



the partition function may be written in a form containing 
a static, i.e. qo dependent, part and a dynamical factor, 
namely, 



Z(0) 



Here 







2ir\k\ J_ 



dq exp[-/3^ A (/3, go )] C(p,q ). 



^ FA (/3,go) 



2|fc 



ql 



lnz(/3,g ) 



(5) 



(6) 



h = h 



FF 



(1) 



Here H and F are hermitian one body operators and the 
product FF mimics an effective separable two body inter- 
action. In the sequel, for the general analysis we always 
assume [H, F] ^ 0, the simplification to the case of com- 
muting operators is straightforward. Further, the coupling 
constant k = — \k\ is taken to be negative (attractive inter- 
action) Jl8| according e.g. to isoscalar modes. The effect of 
repulsive interaction has been studied in The ansatz 
(|l|) defines a minimal microscopic model for a system with 
one collective degree of freedom ]l8| . An extension of the 
methods proposed in this paper to systems with more (in- 
dependent) collective degrees of freedom is then feasible, 
though it may be tedious in detail. 

The partition function of the grand canonical ensemble 
reads 



(2) 



plays the role of an effective static free energy. It is not 
the free energy F((3) — — T In Z(j3) of the total self-bound 
system, but the one of the constituents moving in a mean 
field that is kept fixed at the static collective variable qo. 
For this reason we call ^-" SPA (/3, qo) the "intrinsic free en- 
ergy" from now on, where the index "SPA" already refers 
to the simplest "static path approximation" to (J5|) [see 
Sec. 2.1]. In (|^) there appears the grand canonical parti- 
tion function 

z(P, qo) = Tr exp (-f3(h (qo) - fiA^j 

= Y[(l + exp(-l3(ei(q )-n))) (7) 
i 

belonging to the static part of the Hamiltonian (jl]) in mean 
field approximation 



h (qo) =H + q F. 



(8) 



Z{(5) = Tr exp -/3(H - fiA) = Tr U 



where [3 — 1/T is the inverse temperature (in units with 
fcs = 1) and the chemical potential /i keeps the particle 



We note that our notation has been adapted to the one 
used in (2l||l|] and relates to that of in the following way: 
Ho <-> K, F <-> V, fc < <-> -x < 0, q X <?, z <-> Co, C <-> Co, 

h <-> l. 
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Obviously, /io(<7o) is a sum of only one body operators with 
eigenenergies e;(<?o)- 

All contributions from the dynamical part of the aux- 
iliary path q(r) are contained in the factor C(/3, qo) which 
can be formally written as H 



f(. i.q„) = I V'q expf 



Ilia 



(9) 



with the measure 



(JV-l)/2 



V'q= lim TT -£-r dRe(q r ) dlm(q r ) . (10) 

Ne=hp r=l 1 1 

Here, the thermal expectation value of the evolution op- 
erator U q has to be evaluated with respect to the static 
one-body Hamiltonian (^J) and reads in terms of a time- 
ordered product 



(11 



g/Qo 



f exp [ -- 



dr /ii(r, q r 




(11) 



The "dynamical" Hamiltonian 

h 1 {r,q r )=F{r)8q{r) (12) 

with Sq(r) — q(r) — qo may be understood as the time 
dependent correction to the static mean field Hamiltonian 
^o(<Zo) given in (||). Hence, it is convenient to work in 
an interaction picture based on ho(qo) and define time 
dependent operators as, e.g. 

F(t) = ^ l(, ( qa '> T / h F e ~' lo ( 9o ) T/ ' ?i . (13) 

The partition function (||) together with (0) and (||) is 
still exact, however, written in a way which allows for a 
systematic approximative evaluation. Namely, all quan- 
tum fluctuations, i.e. modes with Matsubara frequencies 
v r 7^ 0, are hidden in (U q ) qa . Thus, the basic idea is to 
successively account for dynamical information in Z((3) 
by expanding the thermal expectation value (l4q} qo around 
its static value in terms of the deviations Sq(r). This way, 
generalizing Eq. (27) of Q| to fourth order in the q r one 
finds: 

inH)° 9 r spA 

^E/ IrQse^e^ (fF(r r )F(r s )) go 



2! h 2 

-y 



3! h? 



hi E 



4! ft 4 



q r q s q t e r r e " s e 1 1 



X (fP(T r )F(T s )F(T t )), 



q r q s qtq u e'^e^V^e 11 ^" 



x (TF{ Tr )F{T s )F{T t )F(r u )) 



8 ft 4 

x (fF(r r )F(r s )) 90 (fF(r t )F(r n )) 90 
0(g r 5 ) (14) 



The symbol / abbreviates summation over all involved 
r, s, . . . ^ and integration of all involved r r , r s , ... from 
to h/3. Because of the r-integrations and the fact that 
(F(r)) qo is r-independent all terms involving such a factor 
— for instance terms linear in q r — vanish immediately 
and are therefore omitted in (|14|). 



2.1 Conventional PSPA - Expansion to second order 

The simplest approximation to -Z(/3), coined the static 
path approximation (SPA), neglects all dynamical contri- 
butions, i.e. one puts hi (r, q r ) = 0. As a result C SPA (/3, go) = 
1. The SPA is the classical limit where all many body 
quantum fluctuations are absent. For lower temperatures 
when quantum properties tend to become important, fluc- 
tuations around the classical limit can be incorporated 
within the conventional version of the perturbed static 
path approximation (PSPA). There, the expansion rtlj) 
is truncated after the second order terms in the q r []6£Q, 
§,H0,0,H|lIl|iJ, which effectively means to describe 
quantum effects on the RPA level. Within the PSPA the 
dynamical factor C(f3,qo) is approximated by |12|] 



C PSPA (P,qo) 



(15) 



J V'q exp ^--^^ {1 + kx(qoAvr)) q r q-,^j 

Here, the FF-response function x(go,cj) defined by 
£(^)<2oH = -x(9o,w) 6q(u) 



(16) 



has to be evaluated at the Matsubara frequencies \v r along 
the imaginary axis. Because of q- r = q* all integrals in 
( |l5| ) are of Gaussian type, and cause no problem as long 

as 

XrW,qo) = l + fcx(?o>r) > for all r >0 . (17) 

For systems with go-regions where the local RPA frequen- 
cies LUv(j3,qo) obtained by the condition l + kx(qo, u v ) = 
are purely imaginary eq.(^) defines a condition for that 
temperature below which the conventional version of the 
PSPA breaks down due to a vanishing Ai(/3, qo). This 
breakdown temperature To where Ai(l/To, qo) = for the 
first time is known from the theory of dissipative tunnel- 
ing as the "crossover temperature" (for many body 
systems sec [[12[) 



max 



fc MI 



2tt 



(18) 



with wJ, nst (<7o) as imaginary local RPA frequencies. Now, 
the condition (|l7|) is guaranteed for all T > Tq and we ob- 
tain for the partition function within conventional PSPA 
[0,|l^] the integral (||) with the dynamical factor C(/3, qo) 
providing improvement over pure SPA taken as 



C 



PSPA 



(A«b)=n 



i 



r>0 



K{f3,qo) 



(19) 
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2.2 Extended PSPA - Expansion to fourth order 

The instabilities that lead to the breakdown of the conven- 
tional PSPA at To are due to the fact that all fluctuations 
q r are treated only up to second order. Since the Gaussians 
in q r have a typical width [see (na)] 



1*1 



\k\T 



2/?A r (/?,g ) V2A r (/?,g ) 



(20) 



first Ai grows as temperature approaches T from above, 
i.e. Xi(/3,qo) — ► 0. Correspondingly, fluctuations in qi- 
direction increase which in turn renders the Gaussian ap- 
proximation for the gi-mode insufficient. In contrast, the 
other amplitudes q2, 93, ... remain still sufficiently small 
for a harmonic approximation to be valid. Therefore in 
the expansion (Jl4|) terms higher than second order must 
be taken into account only for the q\ and g_i mode. This 
is very similar to the procedure used in the framework of 
dissipative tunneling in order to overcome the irregularity 
of the decay rate at To J2^]. Working along these lines a 
consistent approximation to the partition function (|5|) is 
obtained for temperatures near T by taking for C(f3, go) 



C 



cPSPA 



J ^'q exp 



A{(3,q ,q T )\ , (21) 



The relevant effective multidimensional fluctuation poten- 
tial turns out to be 



A(J3, go, q r ) = 2J MA Qo) Qrq- 



(22) 



r>0 



+3c^"(/3, q ) qfq-2 + 3c 3 ((3, q ) q\q 2 
+6c 4 (/3, g ) <??<?-! • 

Here, other terms e.g. of the type gig s g_ s _i or gig_ig s g_ s 
(s > 1) are negligible in the sense of a semiclassical ap- 
proximation as they contribute only in lower order j22j. 
The relevant coefficients containing third and fourth or- 
der F correlations read: 



4 (A qo) 

~\k\/P 

31 h 3 

(/3,go) 

3! h? 
4! ft 4 



1 TV 1 T s W_ 2 T t 



^iiz-it,. ii/_it s ii/ +2 r t 



(23) 

(fF(T r )F(r s )F(T t )) qo 
(24) 

(fF(r r )F(r s )F(r t )) qo 
(25) 

JV+ 1 Tr & W+ 1 T s gi^- 1 T t 1 T u 

((fF(T r )F(T s )F(T t )F(T u )) qo 

- 3(fF(r r )F(r s )) qo (TF(T t )F(r u )) qo 



The integrals symbolize integrations over all involved r 
from to Kf3. Mind that only Matsubara frequencies with 



indices r,s,t,u = ±1,±2 enter these formulas. The ex- 
plicit evaluation of these coefficients is a crucial step in 
order to apply the ePSPA to man y bo dy systems as we 
will demonstrate in d etail in Sec. 4T. In particular, we 
will prove in Sec. 4.1 [see (fi"6|)] that the r integrations of 
( |i~4| ) imply the sum rule r + s + 1 + . . . = for the indices 
involved which considerably simplifies the calculation. 

Let us now turn to the q r integrals in (pl|). Integrals 
with |r| > 1 are still of Gaussian type and can easily be 
carried out as long as A2(A qo) > is sufficiently large. 
While all integrations for |r| > 2 just provide factors 1/A r 
like in conventional PSPA, the integrals over qi and g_2 = 
q 2 require special care. To this end the variables qi and g_2 
are transformed to variables q' 2 — Re q 2 and q 2 — Imq 2 
which cast the relevant part of (S) in the form 



a((3,q ,q' 2 ,q 2 ) = A 2(A<7o) {q' 2 ) 2 
+ 3 (cj(P,q )q 2 _i 
+ X 2 ((3,q )(q^ 2 
+ 3i (t£(p,qo)£ 



(26) 



cf(P,qa)qT, 



q'2 



4- 



The corresponding two Gaussian integrals can be easily 
solved: 



W\ I dq ' 2<lq[ 



2 ex P 



— a{(3,q a ,q' 2 ,q' 2 ') 



(27) 



1 P1fn (£ 9c 3 " (A g Q )4 (A go) 22 
A 2 (A*)) P \\k\ A 2 (A<Zo) 



What is left are the gi-integrals. These are integrals 
over an exponential with an exponent given by 



M(i3,qo)qiq-i+ B(f3,q )qfq 2 _ 1 , 
where the coefficient B(/3,qo) is defined as 

9c3 (P,qo)4 {P,qo) 



B(P,qo) = 6c 4 (/3, g ) 



A2(/3,go) 



(28) 



(29) 



To proceed further it is convenient to introduce polar co- 
ordinates p and 4> instead of gi and g_i = q*, such that 
the phase integral can be done. After an additional sub- 
stitution z = p 2 the gi integrals reduce to 



JL 

n\k\ 



J d<f>J^ dpp «q>(-|fi (Aip 2 + Bp 4 )) (30) 

(3 f°° ( (3 

W\Jo ^ e ^{-Jk\( X ^^ z + B 



(P,qo)z" 



Eventually, together with the contributions from (|27^ 
and ( |30l) the expression ( ^l|) for the new dynamical factor 
reads in "extended PSPA" 



C cPSPA (/3,g ) = n 



1 



dz exp 



r> ^ K([3,q ) 

(3 



(31) 



1*1 



(X 1 (p,q )z + B(P,q Q )z 2 



In case that gi fluctuations are small so that third and 
fourth order terms in gi and g_i are negligible — -B(A go) <SC 
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Ai on) — the factor C oPSPA (/3, q ) reduces to C PSPA (/3, q 
[see (U9h]. For temperatures T < Tq and go for which 
Ai < but | Ai| again sufficiently large, the integral in 
( |3l| ) can be evaluated in stationary phase approximation. 
The result is 



c cpspA ( / 3, (Z o)=n 



1 



\ A r (/3,go) 

27T/3 

\k\B(P,q ) 



(32) 



exp 



d Af(/3,g ) 
\k\4B(/3,q Q ) 



Here, the exponential reveals the change of stability at 
T = Tq as now the contribution to Z((3) of those go 
where Ai < is enhanced compared to their static value 
exph/3^ SPA (/3,g )]. 

The ePSPA smoothly connects the temperature range 
above T with that below Tq. Its precise lower bound of 
validity is determined by two conditions: B(/3,qo) defined 
in (29) and A2(/3, go) [see (|l7|)1 must both be positive and 
sufficiently large for all go- Even if B > also for T < T 
the ePSPA definitely fails at T = To/2 where A2 vanishes. 
Physically, the instability of the PSPA at T = T corre- 
sponds to an instability in gi direction in functional space 
of that classical mean field solution go = g§ for which 
Ai(l/To,g§) = 0. At the temperature where the ePSPA 
breaks down, i.e. where A2 = 0, the classical field at gg 
becomes unstable in g2 direction in functional space. This 
scenario proceeds with decreasing temperature. 



2.3 Treatment of low temperatures T< T 



In the last section we have proposed an extension of the 
conventional PSPA that is applicable down to tempera- 
tures not too far below To. As already mentioned, this ap- 
proximation breaks down near A2 = where the g2-mode 
amplitude becomes large. Naively, one could think to "reg- 
ularize" the g2-mode divergency in the same way as for the 
gi-mode around Tq. However, if both gi and g2 are large 
even higher order coupling terms to q r modes (|r| > 2) in 
the expansion must be taken into account [compare our 
remark after eq. (p2|)j . This proceeds to the divergencies at 
A3 = and so on. Hence, practically an analytical treat- 
ment analogous to the To case is no longer possible for 
temperatures T <C T). In particular, for T — > all modes 
q r must be assumed to be large so that in the worst case 
any kind of semiclassical approximation to Z(f3) fails. 

On the other hand, for T < T but still T > 0, we 
expect that integrals of the type (||) are dominated by 
the contributions around the minima g™ n , i = 1, . . . , M 
of the intrinsic free energy (||) provided /?.F SPA (/3, g^) ex- 
ceeds (3!F SPA ((3, q™l n ) by terms sufficiently larger than of 
order 1 . Under the additional condition that the correction 
factor C ePSPA (/3, go) is a smooth function of go compared 
to the exponential we can apply a saddle point approxi- 



mation to evaluate (||) together with (|3l|): 

1 



M 



ePSPA 



(33) 



t= ^|A|d^A /dg g| cr 

x exp [-/3JF SPA (/?,g™")] C cPSPA (/3, q^) ■ 

The assumption that C ePSPA (/3, go) is smooth compared 
to exp[— /?.T SPA (/3, go)] — variations of the dynamical fac- 
tor are negligible on the scale on which the exponential 
varies — bases on the fact that all approximations of SPA 
type are semiclassical approximations. This condition is 
certainly not fulfilled in regions where the dynamical fac- 
tor diverges. The divergencies, whenever \ r {/3,qo) = 
(r > 1) only appear in those go-regions where the RPA fre- 
quencies are imaginary and of sufficiently large amount, 
thus, corresponding to saddle points of the full free en- 
ergy. Therefore we may assume that they need not be con- 
sidered for low temperatures (but still sufficiently above 
T = 0). Accordingly, the following interpolation between 
the ePSPA and the saddle point approximation ( |33| ) seems 
to be natural: For T <C T) we restrict the go integration 
in (H) with (|3l| ) to those regions where Ai > meaning 
that all X r > 0, i.e., 

C LTA (/3,g ) = 0[Xi(P,qo)} C cPSPA (/3, g ) (34) 

where #(•) denotes the step function. This interpolation 
is called "low temperature approximation" (LTA) hence- 
forth. It requires to solve Xi(f3, g§) = for given (3 > (3q = 
1/To to determine g§(/3). As a continuous function of (3 
starting at gg(/?o), go(/3) describes the broadening of the 
instability region around go(A)) where Ai < |23|j . The 
ePSPA extends to temperatures somewhat below To, while 
the LTA is valid if with decreasing temperature T < To 
the contribution of the instability region around go (A)) to 
the go integral becomes negligible. Therefore both approx- 
imations match in a narrow temperature range below Tq. 
Further, by formally extending the LTA to T > To one 
recovers the ePSPA since then Ai > for all go, while for 
T ^ Tq the saddle point result ([33]) is regained. As al- 
ready mentioned, in general there is a lower bound for the 
LTA. Contributions of mean fields around g§ can only be 
neglected as long as the exponential enhancement in (|3^ ) 
due to the dynamical instability is still compensated for 
by the difference between F SPA {f3, g§) and F SPA {f3, g™«). 
For very low temperatures this can no longer be taken for 
granted. 

We want to remark that for specific systems further 
approximations may exist in the region T <C Tq, possibly 
better adapted to the problem in question. This may be 
particularly true for very low temperatures T — * if new 
phenomena, sometimes associated with the appearance of 
new symmetries like in the case of Goldstone modes, re- 
quire special care (cf. the cases of a one-dimensional dou- 
ble well and the Lipkin-Meshkov-Glick model below in 
Sects. [|and^|). As long as one is interested in approxima- 
tion schemes applicable to general situations without re- 
ferring to individual system properties, however, the LTA 
presented above seems to us the only consistent one. 
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3 Application I: One dimensional potential 

In order to illustrate the utility of the the above approxi- 
mations we turn to exactly solvable models and first apply 
the formalism to the one-dimensional case. 

Consider a particle of mass M in a one dimensional 
potential V(q). The imaginary time path integral repre- 
sentation is given by 

Z{p) = J dq' j>Vq(r) exp (-S E [q(r),q(r)}/h) (35) 



with the Euclidean action 

cfifi 



S^q(r),q(r)} = 



M 



q 2 (r) + V(q(r)) 



(36) 



In ( pq ) we first have to sum over all paths q(r) starting 
and ending at a given end-point q(0) = q{h(3) = q' and 
second sum up contributions from all these end-points q' . 

In order to be able to apply the ideas of Sec. the 
potential must be expanded up to fourth order in the 
Fourier coefficients q r (\r\ > 0) introduced in (||). After 
carrying out the r integration we obtain for the exponent: 

- S E [q r ]/h = -PV(qo) {Mvl + V"(q )) \q r \ 2 



3! 
4! 



r>0 

S r +s+t,o q r q s qt (37) 

r,s,t^0 

3r+s+t+u,o qrqsqtq-u ■ 



For the case of one dimensional potentials the sum rule 
r + s + t + u = simply stems from this r integration. Here 
V(qo) plays the same role as F SPA {f3, q ) in (|). As before 
the coefficients q\ and <?_i become large at the crossover 
temperature 



To 



h_ 



-V"(q ) 



AI 



(38) 



and therefore have to be taken into account up to fourth 
order. Then, with \ r (f3,q ) = Mv 2 + V"(q Q ), c+(f3,q Q ) = 
c-(l3,q Q ) = ^ 3 )(<z )/3!, and c 4 (/3,9o) = V^(q )/4l we 
gain 



AI 



2Trh 2 f3 



dq exp[-(3V(q )] C(f3,q ) . (39) 



- F(P) 

A -HQ 




** - \ 

\ 
1 




exact 




SPA 




PSPA 




ePSPA 




LTA 



2 4 ( 

Fig. 1. of the double well 

approximations and A = 0.9. 

and for the extended PSPA 

£tePSPA 



8 10 

2h computed for various 



dz exp < — j3 



2 r + V"(q ) 
Mv\ + V"(q )) z 



(41) 



Let us now turn to the specific case of a particle of 
mass M = 1/2 {%= 1) moving in the quartic potential 



V(q) = ~Aq 2 + {1- A) q 4 + 



A 2 



4(1 -A) 



with < A < 1 
(42) 



This potential coincides with a quartic oscillator for A = 
and develops a barrier of height V(0) = A 2 / (4(1 — A)) for 
positive A. Potentials of the form ( fi"2"|) have often been 
used to qualitatively understand characteristic properties 
of high dimensional systems as e.g. phase transitions. Here 
we find that the shape of V(q) describes the main features 
encoded in the intrinsic free energy of the Lipkin-Meshkov- 
Glick model which we will consider in detail below (see 



Sec. 4.2). There, the appearance of a phase transition is 
equivalent to the development of a barrier for A > in 
V(q). Further, from pi] ) and ( [42] ) the inverse crossover 
temperature can simply be derived as /3q — 1/Tq = 2ir/h- 
y/M/(2A) = n/(hVA). In fig. g the free energy 



1 



H3)=--ln£(/?) (43) 
is depicted within various approximations for A = 0.9. 



In case of the pure SPA we have C SPA (/3, qo) 
of the conventional PSPA 



c pspA (/3, 90 )=n 



Muz 



r>0 



Mv?+V»(q ) 



1, in case 



(40) 



4 Application II: The Lipkin-Meshkov-Glick 
model 

The approximations derived in the preceeding sections 
are designed to describe low temperature thermodynamic 
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properties of interacting many body systems. In contrast 
to simple one-dimensional cases, here, an additional chal- 
lenge arises through the calculation of the expansion co- 
efficients in (|l|) and 



4.1 Expansion coefficients for many body systems 

For the Hamiltonian ([!]) the coefficients ct and C4 are 
calculated by evaluating the r-integrals in (p3|)-(p5|). The 
calculation is performed completely parallel to the one 
presented in |7|. Hence, we restrict ourselves here to one 
contribution to the fourth order coefficient to exemplify 
the general strategy. We also will prove the sum rule r + 
s + t + u = exploited in ( |22| ) ( pq ) here. To this end we 
consider general indices r, s, t, at first. One has to 
rewrite the operators F(r) in terms of creation and an- 
nihilation operators F(t) = J^jy ^fej^i( T )% ( r ) an d then 
calculate the time ordered averages using the finite tem- 
perature Wick theorem Evaluating all possible con- 
tractions the unperturbed temperature Green functions 
enter, i.e. 



(Ta^r^aliTs))^ = S jk g^(r r - r s ) (44) 

1 J2, g-iWA-(T r -T 3 ) 



,(0), 



(3 ihiu K - e k (q ) 



where u>k = {2K + 1)tt lh(3 (^|) and efc(go) corresponds 
to the eigenenergies of (g) via efc(go) = efc(<7o) — M- The 
integrands of all multiple r integrals therefore turn into 
factors like e.g. 



exp 



1 — (r + I 



hp 



0)t 7 



(45) 



After carrying out all r integrals we obtain the conditions 
r + I - O = 0, s - I + K = 0, t- K + M = and 
u + O — M = 0. Putting these conditions together we 
easily obtain the requirement 



r+s+t+u=0 



The infinite sum over I is calculated by exploiting the fre- 
quency summation technique One uses the replace- 
ment icjj — > z and the residue theorem for contour inte- 
gration in the complex z-plane together with the obser- 
vation that the function —hf3/ (exp(h(3z) + 1) has poles at 
iwj with residues 1. After a deformation of the integration 
contour such that it encycles the poles of ( [47j ) we obtain 
the following final form: 



\k 



~~M E F io(Qo)Fki(qo)F m k(qo)F om (qo) 



(48) 



i.k.m.o 



1 



1 



1 



+ n(e ) 
+ n(e k ) 
+ n(e m ) 



e lo + \hv r e ik - \hv s e im - \hv s+t 
1 1 1 



£oi — ihi/ r e k — ifev+s e om — ifitv+s+t 
1 1 1 



tki + ifiv s eko + iftvr+s e km - ifiut 
1 1 1 



Here n(e) = (l+exp(/3(e— fx))) 1 are the Fermi occupation 
numbers and Cik{qo) are the energy differences between the 
states \i(qo)} and \k(qo)) of the static Hamiltonian ho(qo). 
This way, by solvingthe static one body Schrodinger equa- 
tion belonging to (gj) and evaluating sums as in ( [48| ) all 
needed coefficients are known and the corresponding ap- 
proximations to Z((3) can be applied. 

We note that the factor in brackets in (|48| ) also ap- 
pears in |tJ where strength functions are calculated in 
the conventional PSPA formalism. There, time ordered ex- 
pectation values of products of one body operators D = 



kj ^kj a Jfe a j 



this end U q [see 



like {TU q D^ (t)D{Q)) qa are evaluated. To 
HI) and (14)] have to be expanded up 



to second order in the q r . As a consequence, fourth or- 
der terms in a} k {T)aj(r) come into play even though the 
approximation to lA q is still of Gaussian type. 



(46) 4.2 The model 



Moreover, we have only one summation index / = —00 ... 00 
left besides r, s, t 7^ which will actually be further re- 
stricted by the requirement r, s, t, u = ±1,±2 in ( p3| ) - 
( p5h . One can show after some straightforward but lengthy 
algebra that one typical contribution to the fourth order 
coefficient C4(/3, go) reads: 



4! 



There are only few many body systems which allow for an 
exact evaluation. One is the Lipkin-Meshkov-Glick model 
[ p~7| which has been used to test the results of conven- 
tional PSPA for years (see e.g. ||,@,§|). In the last part of 
the paper we want to apply the approximations derived in 
Sec. U to this model. The Lipkin-Meshkov-Glick Hamilto- 
nian reads 



E 



Fio(qo)Fki{qo)F m k(qo)F om (q ) 



H = 2eJ z + 2Ul 



(49) 



1 



i = -c 



Kusi 



ei(go) ih(u!i 
1 



e Q (go) 
1 



(47) 



\h(ijji 



ek(qo) - v„ +t ) 

-esp 

are equivalent to the ov of |7|. 



Mo) 



2 Mind that our uk correspond to the v k of and our 



and has the structure of a Hamiltonian with two body 
interaction with H = 2eJ z and F = 2J X . In ( |4g| ) the op- 
erators J x and J z obey angular momentum commutation 
relations and k = —\k\ is a negative coupling constant 
describing an attractive interaction. 

Let us briefly recall ||[7],|j| the most important features 
of the PSPA applied to the Hamiltonian (49). The static 
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part (^|) of the one body Hamiltonian is given by 

h (qo) = 2eJ z + 2q J x (50) 
and its eigenvalues are 5-fold degenerated 

e 2 M = e 2 +<Zo • (51) 

Energy differences tik{qo) can therefore only have the three 
different values 0, ±2e(q r o). For the grand canonical parti- 
tion function belonging to ( [sof ) one easily obtains at a 
given collective coordinate qo [see ([?])] 



2 cosh 



2;/ 



The PSPA correction factor ( [19|) can be written as 
sinh(/3e(g )) Hf3w(/3,q )/2 



C 



PSPA 



f3e(q ) emh(h0m(ft, q )/2) 
where the RPA frequencies read 



(52) 



(53) 



/ Hw(P,q ) 



= e 2 (qo) - k 



e (?o) 



tanh 



/3e(9o) 



with the dimensionless coupling parameter (| 



(54) 



(55) 



For k > 1 the equilibrium value of q$ undergoes a phase 
transition at a critical temperature T c which is deter- 
mined from the condition that one can find a go with 
zero RPA frequency from (^4j) . This phase transition man- 
ifests itself in the development of a barrier in the intrinsic 
free energy for T < T c corresponding to imaginary RPA 
frequencies. Thus, the crossover temperature To where 
the factor (|5^) diverges is determined from the condi- 
tion (hf3m((3,qo)/2) 2 — — 7r 2 . We note that one always 
has To < T c . For k < 1 no phase transition occurs and the 
RPA frequencies are always real. 



4.3 Free energy 

For the Lipkin-Meshkov-Glick model the free energy has 
been calculated in j7j] within the conventional PSPA. Here, 
we want to demonstrate the results of the extensions pro- 
posed above using the same set of parameters as in ||, 
namely, 



e = 5 MeV and g = 10 



(56) 



and various values for n. 

In figs. we present plots of the free energy ( (43| ) 
as a function of the inverse temperature f3. Fig. @ shows 

3 Our k is defined as in |7J and should not be mixed up with 
the k = 1/k often used in nuclear physics Jlit - 
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Fig. 2. T{f3) of the Lipkin-Meshkov-Glick model ([49|) within 
various approximations for the parameters ( |5q ) and k — 1.313. 
The phase transition occurs at f3 c — 0.4 MeV - and the inverse 
crossover temperature is fta = 1.1392 MeV -1 . 



the global behavior for k ~ 1.313 (the value used in M). 
Figs. ||a-c illustrate a blow up of the region around T for 
three different coupling strengths. The exact result can be 
obtained by a numerical diagonalization of the Hamilto- 
nian (M9I) in a basis of eigenstates of J z or J x as explained 
in [0,0)0 • As expected, pure SPA deviates strongly from 
the exact results with decreasing temperature. For n < 1 
no phase transition occurs and the difference between con- 
ventional and extended PSPA is very small [see fig. 0a]. 
The inclusion of quantum effects on the RPA level in the 
conventional PSPA (|9|) improves the results a lot as long 
as (3 is not too close to the breakdown value /?o = 1/Tq 
[see fig. |[. Here the conventional PSPA fails as the correc- 
tion factor diverges for k > 1. Instead, the extension ( |3l| ) 
smoothly passes the crossover region and agrees very well 
with exact results even for inverse temperatures slightly 
above (3q. The ePSPA definitely br eaks down near (3 = 2(3q 
where A2 = 0. The LTA [see Sec. 2.3 smoothly matches 
with the ePSPA somewhat below (3 and is able to cover 
even the range /3 3> /3o with astonishing accuracy. 



One can also see in fig. |3Jb that the matching between 
ePSPA and LTA even though it appears to be continu- 
ous in temperature, may lead to discontinuities in first 
and higher order derivatives. The latter ones produce un- 
physical results for the internal energy and specific heat 
(see below) around the matching point close to To- This 
kind of behavior is typical for the matching between dif- 
ferent semiclassical approximations each one designed for 
different ranges in parameter space. It can be overcome 
in principle by invoking uniform semiclassical approxima- 
tions. Apart from that the ePSPA together with the LTA 
provide a practicable and very efficient way to evaluate 
the free energy over a broad temperature range. 
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Fig. 3. Blow up of the crossover region for T{(3) of the Lipkin-Meshkov-Glick model (|49|) within various approximations. 
Parameters are as in (^6|). For k — 0.8 no phase transition occurs: To = 0. For « = 1.313 the inverse crossover temperature is 
f3 = 1.1392 MeV" 1 whereas for k = 1.8 the corresponding value is (3 = 0.7433 MeV -1 . 
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Fig. 4. £(J3) of the Lipkin-Meshkov-Glick model @ for the 
parameters of ( |5rj| ) and k = 1.313. The phase transition occurs 
at f3c = 0.4 MeV -1 and the inverse crossover temperature is 
0o — 1.1392 MeV -1 . 
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Fig. 5. C(/3) of the Lipkin-Meshkov-Glick model @ for the 
parameters of ( |5rj| ) and k = 1.313. The phase transition occurs 
at f3 c = 0.4 MeV -1 and the inverse crossover temperature is 
j3 = 1.1392 MeV -1 . 



4.4 Internal energy and specific heat 



4.5 Level densities 



The internal energy (comp. g) 



d 



(57) 



is shown in fig. [|. In the crossover region (3 w /3o the ex- 
tended PSPA again is an improvement over conventional 
PSPA but tends to fail for larger (3. The low temperature 
approximation of Sec. 2.3 delivers reasonable values for 
[3 > Po. As already mentioned above, it shows unphysical 
discontinuities close to /3q — 1/Tq. 

The specific heat of the system (comp. H) 



C((3) 



(58) 



is seen in fig. |5[ Pure SPA turns out to provide a reason- 
able global approximation of the specific heat over a broad 
range of temperatures, whereas the conventional and ex- 
tended PSPA lead to very good results at higher temper- 
atures, smaller f3, but deviate already for (3 < (3 . The 
LTA of Sec. 2.3 again shows a discontinuity at f3 = f3o and 
supplies a reasonable estimate for (3 > (3q. 



As already mentioned in the introduction the develop- 
ment of the conventional PSPA was motivated to some 
extent by the need to calculate nuclear level densities as 
a function of the excitation energy £*. Given the parti- 
tion function this can be achieved by an inverse Laplace 
transform. Within a saddle point approximation (Darwin- 
Fowler method) one finds the inverse thermal temperature 
(3* from £ * = £{(3*) — £q and arrives at the following for- 
mula for the level density p5[: 



Pin 



i 



exp[«S(/3*) - (3*£ ] 



The entropy is given by <S(/3) = (3 [£{(3) — J-(0)] and 



V 



06 



2 lnZ(/3) 



0=13' 



1 



03* 



■C((3*) 



(59) 



(60) 



Results are depicted in fig. q Here, high excitation ener- 
gies correspond to high temperatures. For this reason ap- 
proximations of SPA type are expected to be exact at large 
£* as confirmed in fig. pi For low excitations the situation 
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is different. The pure SPA is not able to deliver results for 
£* < 4 MeV. The reason is simply that the SPA result for 
the internal energy is about 4 MeV larger than the exact 
result at small temperatures [see fig. || and therefore the 
condition £(/3*) - £ = £* < 4 MeV cannot be fulfilled 
within SPA. The conventional and extended PSPA results 
agree with the exact ones much better than those of pure 
SPA. Essential deviations only occur for £ * < 2 MeV and 
the extended version is a little better than the conven- 
tional one [see fig. m. 

The main problem of the conventional and extended 
PSPA with respect to the level density can be traced back 
to the quantity V of ( |60| ) that is proportional to the spe- 
cific heat C(/3*) and enters in the denominator of (p9|). 
Divergencies in the specific heat as seen for PSPA and 
ePSPA in fig. |5| therefore have a tremendous impact on the 
results for the level density. Therefore both types of PSPA 
deliver level densities that are too small which becomes 
significant at small excitations corresponding to (3 > [3q. 
The following strategy seems to be reasonable: In order to 
calculate the level density using (|5^) we exploit an approx- 
imation for the specific heat that is globally reasonable and 
combine it with the best approximation at hand for the 
internal and free energy. That means to use pure SPA (see 
fig. |) for the specific heat and the ePSPA/LTA for the 
remaining quantities. As shown in fig. [?] this procedure 
improves the conventional and extended PSPA a great 
deal at low excitations. Remarkably it even describes the 
bending up of the level density qualitatively correct. 



5 Conclusion 

We have applied the path integral approach to approx- 
imately evaluate the partition function of a finite inter- 
acting many body system in the low temperature regime. 
This requires to extend the conventional PSPA down to 
temperatures below the crossover temperature To. At these 
temperatures the simple mean field solution becomes un- 
stable and large quantum fluctuations arise. The crucial 
step is to go beyond the Gaussian approximation on which 



the PSPA relies and take into account anharmonicities of 
certain critical vibrational modes. This procedure stabi- 
lizes the semiclassical type of approximation to the path 
integral of the partition function and holds true down to 
temperatures T rs T)/2. There, an instability in another 
direction in functional space shows up, which successively 
proceeds to occur at all T r = T)/r (r > 1). In order to 
treat the temperature range T <^ T), but still sufficiently 
above T = 0, we proposed an approximation where only 
mean fields with small quantum fluctuations around the 
minima of the static free energy are taken to contribute. 

In this way, we could study thermodynamic properties 
of the archetypical Lipkin-Meshkov-Glick model even far 
below the crossover temperature. In particular, for the free 
and internal energies results were gained, that agree very 
well to the exact ones. Combining our extensions to the 
PSPA with pure SPA we obtained very accurate results 
for the level density even at low excitation energies. 

Improvements beyond the PSPA for T sa T and T < 
T) are desirable in order to get reasonable approximate re- 
sults for many body systems in the low temperature range 
where usual Monte Carlo techniques become very expen- 
sive. This is particularly true as often the critical tem- 
perature T c where a phase transition takes place and the 
crossover temperature Tq where the PSPA breaks down 
are of the same order of magnitude. For instance, for finite 
nuclei one finds for the nucleus 164 Er that T c « 0.45 MeV, 
while PSPA is reliable only for T > 0.25 MeV Q. In |oj 
the range of applicability of PSPA is given by T > 0.2 MeV 
for the nuclei 104 Pd and 114 Sn. 

In general, for these or even lower temperatures (of 
the order of 10 _1 MeV and below for nuclei) the concept of 
temperature is at least doubtful for small isolated systems. 
This is different for systems on a mesoscopic scale. There 
the coupling to a macroscopic heat bath is always present 
such that temperature can be fixed from outside. Further, 
the number of constituents is typically much larger as e.g. 
for finite nuclei. Accordingly, as addressed above, the con- 
ventional PSPA has been applied to small superconduct- 
ing and superfluid systems. For superconducting particles 
of nanometer scale one roughly finds Tq ~ T c /2 |H| so 
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that an extension of the PSPA is clearly relevant. Work 
in this direction is in progress. 
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